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^ : ABSTRACT 



O ! We report the first high-resolution X-ray spectra of two high-redshift quasars, 

S5 0836-1-710 and PKS 2149-306, obtained with the Chandra High Energy Trans- 
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^ \ mission Grating Spectrometer (HETGS). The primary goal of this observation is 

to use the high spectral resolving power of the HETGS to detect X-ray absorption 
produced by a hot intergalactic medium. The continuum of both quasars can be 
fitted by absorbed power laws. The power law photon index (F) of S5 0836+710 
is consistent with that found in previous observations with ROSAT and ASCA, 
^ \ while the power law of PKS 2149-306 is harder than found in a previous obser- 

vation by ASCA. Excess continuum absorption above the Galactic value is found 
in S5 0836+710, as evidenced in ASCA and ROSAT observations. No significant 
emission or absorption feature is detected in either source at ±3(T level. Based 
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H ' on the detection limits we constrain the properties of possible emitters and ab- 

ed ' 

sorbers. The upper limit of the equivalent width of Fe K emission lines could be 
as low as ~ 10 eV . Absorbers with a column density higher than 8 x 10^^ cm~^ 
for O VIII or 5 X 10^^ cm~^ for Si XIV would have been detected. We propose 
a method to constrain the cosmological parameters (namely VLq and Vtb) via the 
X-ray forest theory, but current data do not give significant constraints. How- 
ever, we estimate that observing ~ 7 bright quasars should give at least one 
O VIII and more O VII absorption lines at 95% possibility. We also find that 
combined with the constraints from the distortion of the CMB spectrum, the 
X-ray Gunn-Peterson test can marginally constrain a uniform, enriched IGM. 
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Subject headings: intergalactic medium — quasars: absorption lines — quasars: 
individual (S5 0836+710, PKS 2149-306) — cosmology: observations — X-rays: 
galaxies 

1. Introduction 

Cosmological numerical simulations suggest that a large fraction of the baryonic matter 
should reside in intergalactic space in the form of a hot/warm intergalactic medium (IGM) 
(see, e.g., Fukugita, Hogan & Peebles 1998; Cen & Ostriker 1999; Dave et al. 2000). This 
hot/warm IGM contains at least 30% — 40% of the baryons predicted by the standard theory 
of big bang nucleosynthesis (BBN) at temperatures between 10^ to 10^ K. Observationally 
little is known about this IGM. Recent observations of O VI absorption lines (see, e.g., 
Tripp, Savage & Jenkins 2000; Tripp & Savage 2000; Oegerle et al. 2000) indicate that these 
O Vl-absorbers could be a significant "baryon reservoir" at low redshift. 

Detecting the X-ray absorption features in the X-ray spectrum of a distant quasar was 
first proposed by Shapiro & Bahcall (1980). They suggested that an intervening uniform, 
hot IGM would introduce an "absorption trough" in the X-ray continuum spectrum of a 
distant quasar via the "X-ray Gunn- Peterson test" . Absence of such an "absorption trough" 
would put constraints on the temperature, density and ionization states of the diffuse IGM. 
Aldcroft ct al. (1994) extended this idea by including photoionization from the X-ray back- 
ground radiation. By using numerical simulations Hellstcn, Gncdin & Miralda-Escude (1998) 
explored the idea that a clumpy, highly-ionized IGM concentrated in the gravitational poten- 
tial caused by dark matter might introduce resonant absorption lines in the X-ray spectrum 
of a background quasar and denoted them as the "X-ray Forest", in analogy to the Lja 
forest in the optical/UV bands. Similar ideas were discussed by Perna & Loeb (1998) and 
Fang & Ganizares (2000) (Hereafter FG). 

A rough estimation shows that even for the largest mass concentrated objects, such as 
clusters of galaxies, the maximum column densities of the absorption ions (such as O VIII 
and Si XIV) are between 10^^ - 10^^ cm~^ and the maximum equivalent width (EW) is far 
below 10 eV (see, e.g., Sarazin 1988; David 2000). All previous observations have been 
limited by the spectral resolving power. For instance, ROSAT and ASCA could not detect 
any spectral feature with an EW smaller than 10 eV. However, the launch of the Chandra 
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X-ray Observatory^ and XMM-Newton make it possible to study the X-ray spectrum of a 
distant quasar with high-rcsohition spectroscopy. The High Energy Transmission Grating 
Spectrometer (HETGS, see Markert et al. 1994; Canizares et al. 2000) on board Chandra 
provides excellent spectral resolution with AE ^ 1 eV around 1 keV. Given such a energy 
resolution, it is possible to study the resonant absorption lines or edges produced by the hot 
IGM in the X-ray spectrum of a distant quasar (Canizares & Fang 1998). 

In this paper we present observations of two high-redshift quasars: S5 0836+710 {z = 
2.17) and PKS 2149-306 {z = 2.34). At redshift higher than 2, these two sources are the 
brightest (S5 0836+710) and the second (PKS 2149-306) brightest known quasars in X-ray 
band'^. Given such a high redshift, the light from the source traverses approximately 80% of 
the Universe, which greatly increases the possibility of being attenuated by the intervening 
systems. Our primary goal is to find if the predicted "X-ray forest" is detectable with the 
Chandra HETGS and if not, what are the hmits of possible absorbers. 

This paper is organized as follows: we describe the observations in §2 and the data 
reduction and analysis in §3. In §4 we discuss the constraints on the IGM, the last section is 
a summary. Throughout the paper unless specifically mentioned, we use a Hubble constant 
of Ho = 70/170 km s-'^Mpc'^. 

2. Chandra Observation 

S5 0836+710 and PKS 2149-306 were observed with the High Energy Transmission 
Grating Spectrometer on board Chandra. Detailed target information is hsted in Table 1. 
The observation of PKS 2149-306 was broken into two parts, with approximately 18 hours 
between them. The HETGS produces a zero order image at the aim point of the focal plane 

detector, the ACIS-S array, with higher order spectra dispersed to either side (for ACIS-S, see 
Garmirc et al. 2000). The moderate energy resolution of the ACIS-S is used to separate the 
overlapping orders of the dispersed spectrum. The HETGS consists of two different grating 
assemblies, the High Energy Grating (HEG) and the Medium Energy Grating (MEG), and 
provides nearly constant spectral resolution (AA = 0.012A for HEG and AA = 0.023A for 
MEG) through the entire bandpass (HEG: 0.8-10 keV, MEG: 0.4-8 keV). 

The ACIS-S array consists of 6 chips (S0-S5); however, due to a temporary malfunction 
of a portion of the on-board readout electronics (the FEP problem, see Chandra Proposers' 



^See http:/ /chandra. harvard.edu/ 

^see the data archive under http://heasarc.gsfc.nasa.gov/ 
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Observatory Guide, or POG), only five chips could be used (S1-S5). Due to the radiation- 
induced charge transfer inefficiency (CTI) the energy resolution of the four front-illuminated 
chips has become a function of the chips' row number, with a better resolution near the 
read-out nodes. Therefore, we moved the nominal aim point closer to the read-out nodes by 
about 3.3 mm for S5 0836-1-710 and 2.4 mm for PKS 2149-306 to get better ACIS energy 
resolution. The resulting ACIS resolution is adequate to allow unambiguous separation of 
the spectral orders. 

Data were analyzed with the standard pipeline for the Chandra HETGS provided by 
the Chandra X-ray Center (CXC)^. We use a combination of Chandra Interactive Analysis 
of Observations (ClAO) VI. 1 and custom routines in IDL. The standard screening criteria 
were apphed to the data. We selected photon events with ASCA grades 0, 2, 3, 4, 6 and 
excluded those with energies above 10 keV. We also removed hot columns and bad pixels 
in each CCD chip. Figure 1 displays the zeroth-order images of both targets in a 50" x 50" 
region. Wc have smoothed the images in Figure 1 with a Gaussian (a = 1.2") for visual 
clarity. Both unsmoothed zeroth-order images are consistent with point sources and can be 
well fitted by a single Gaussian with a ^ 0.38" for S5 0836+710 and a ^ 0.35" for PKS 
2149-306. The light curves derived from the zeroth orders of both sources were constant to 
within statistics. 

S5 0836-1-710 was reported as having a one-side, superluminal radio jet (Krichbaum et 
al. 1990; Lobanov et al. 1998) at a position angle (PA) of —146°. The jet has a length of at 
least 1.5", which is at the spatial resolving limit of the ACIS-S. The head of the jet is located 
between 1 — 2" from center with an opening angle of ~ 2°. We examine this region in the 
zeroth order image and find that the Scr upper limit of the flux from this region is < 0.04% of 
the total flux from the source. Compared with 0.4% and 6% in the jets of 3C273 (Marshall et 
al. 2000) and PKS 0637-752 (Schwartz et al. 2000), respectively, we conclude that no X-ray 
jet was found during this observation. The observation of PKS 2149-306 was separated into 
two parts and the zeroth order images show that the aim points of the two observations 
were shifted by approximately 2 pixels. Wc analyzed the two observations separately and 
summed them at level 2 of the CXC pipeline to get a combined spectrum. 



''See http://asc.harvard.edu/ 
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3. Data Analysis 

3.1. Continuum Fitting 

To fit both MEG and HEG spectra, we ignored the 1 - 2 keV portion of the MEG 
spectrum due to unmodelled systematic calibration errors of ~ 10% — 15% and fit the 
remaining part simultaneously with the HEG spectrum. We present the fitted HEG spectra 
in figure 2. The spectral model is a single absorbed power-law where the absorption cross 
section and metal abundances are from Morrison & McCammon (1983). Table 2 gives the 
best-fit parameters. All errors are quoted at 90% confidence. For PKS 2149-306, since fits 
with Nh as a free parameter do not converge, we fix Nh at the Galactic value. 

S5 0836+710 was previously observed by both ROSAT and ASCA (Brunner et al. 
1994; Cappi et al. 1997). The ASCA spectra gave a photon index (F) of 1.45+[!:J^^ or 
1.32]'^q'q3 between 2 and 10 keV, with free Nh or Nh fixed at the Galactic value, re- 
spectively. Two ROSAT observations showed a steeper spectrum (the photon indices are 
1.511q;i2 and 1.551q;22 between 0.1 and 2 keV, respectively). The photon index in this 
observation (F = 1.388 ± 0.012) is consistent with that of ASCA but is flatter than the 
ROSAT result. We flnd excess continuum absorption (above the Galactic value, at A.Nh ~ 
4 X 10^° cm"'^). The Galactic hydrogen column density is 2.83 x 10^*^ cm~^ with an la error of 
1 X 10^^ cm~^ (Murphy et al. 1996). Excess absorption was also found in ASCA observation 
at ^Nh ~ 8.5 X 10^*^ cm^^ (Cappi et al. 1997). We see a decreasing of excess absorption of 
~ 4.5 X 10^° cm^^ in a timescale of ~ 1.5 yr in quasar frame. It is still unclear the causes 
of the excess absorption (instrinsic to the quasar or due to the intervening systems) and its 
variation (see discussions in, e.g., Cappi et al. 1997; Reeves & Turner 2000). The fiux of S5 
0836-1-710 decreased by a factor of approximately 2 between the two ROSAT observations, 
and the ASCA flux (~ 1.4 x 10~^^ ergs cm~'^s'~^ between 2 and 10 keV) is consistent with 
the lower ROSAT flux. These observations show a flux approximately twice that seen with 
ASCA, which is consistent with the higher ROSAT flux. Given the overlap of the Chandra. 
HETG and ASCA bandpass and the small instrumental calibration uncertainties of Chandra, 
we conclude that the different flux levels are intrinsic to the quasar. 

PKS 2149-306 was also observed by ROSAT and ASCA (Siebert et al. 1996; Cappi et al. 
1997). The fiux measured using ASCA between 2 and 10 keV is 9.94 x 10~^^ ergs cm~'^s~^ 
and is about 30% higher than our result. The photon indices measured using ASCA (1.55^*^0 02 
for detector SIS and 1.54+°;°^ for detector GIS) are significantly higher than our value (F = 
1.255 ± 0.020), which indicates spectral variability in this source. This source was also 
observed during the ROSAT All-Sky Survey, but the small number of counts precluded 
spectral fitting; a rough comparison with ASCA also implies spectral variability (Siebert et 
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al. 1996). 

Optically selected quasars can be divided into two categories according to their radio 
loudness: radio- loud quasar (RLQ) or radio-quiet quasar (RQQ) (see, e.g., Avni & Tanan- 
baum 1986). Both of the sources we observed are RLQs. Previous observations by Einstein, 
ROSAT and ASCA showed that the two types of quasars have significantly different X- 

ray properties where typically RLQs have flatter power-law indices (F 1.6) than RQQs 
(r ~ 1.9) (sec, e.g.. Reeves & Turner 2000). Our observation are consistent with this finding, 
although the spectra we observed are significantly harder than the average. 

3.2. Line Analysis 

We are interested in any possible hne features in the quasar spectra. To identify hne 

features, we need a good measurement of the continuum. We use the above absorbed power 
law models for HEG spectra. For MEG spectra, we use the same models and fit the residuals 
with a 5-order polynomial to account for the remaining uncertainty in the calibration. This 
fitted continuum serves as a reference point for measuring line features. Since the MEG 
spectrum spans from 2A to 26A, the 5-order polynomial will affect spectral features larger 
than 5A but will preserve narrow hne features. The MEG spectra for both sources are shown 
in figure 3. The bottom panel of each plot gives x, and the two dotted lines within the x 
plot correspond to ±3(7 level where a is the square root of the observed counts in each bin. 

If the model fits the observation well, the counts in each bin should follow the Poisson 
statistic, with a mean count from the fitted spectrum. To verify this, we calculate the Poisson 
probability Px{> k) for each bin, where x is the expect mean number and k is the observed 
value. Pj;(> k) should uniformly distribute between and 1, given that the model fits the 
observed spectrum. In figure 4 we show the cumulative distributions of Pxi^ k) for both 
quasars. We also show the cumulative distribution of the uniform distribution as the dotted 
line. We run a Kolmogorov-Smirnov (KS) test to test the null hypothesis that the both 
Px{^ k) are uniformly distributed. The KS test shows that the null hypothesis is accepted 
at a significance level of 0.59 for S5 0836+710 and 0.32 for PKS 2149-306. 

A line feature should at least have a signal-to-noise ratio at or above ±3o" level to be 
identified. Figure 3 shows that most of the bins are well distributed within the ±3(7 level. 
The two exceptions occur at around 22.6A for S5 0836±710, where x = 3.52, and at around 
17.76A for PKS 2149-306, where x = 3.51. Assuming a Poisson distribution for each bin, 
we ran 1,000 Monte-Carlo simulations for both spectra and find that the probabihty for 
observing one bin with x > 3.52 in one observation of S5 0836±710 is 63.5%, and observing 
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one bin with x > 3.51 in one observation of PKS 2149-306 is 48.9%. So we conclude that no 
hne feature is identified. 



3.2.1. Emission Lines 

Iron emission hues have been detected in several high-redshift quasars with ASCA (see 
Reeves & Turner 2000 and references therein). However, no previous observation discovered 
any emission hues in these two sources, with the exception of Yaqoob et al. (1999) (Y99 

hereafter) for PKS 2149-306, which wc will discuss later. In ASCA observations the upper 
limit on the equivalent width of an Fe Ka emission line (~ 6.4 keV) is 110 eV for S5 
0836+710 or 85 eV for PKS 2149-306, respectively (Cappi et al. 1997; line measured in the 
observer's frame, with line width cr = for an unresolved line). 

The Fe K emission lines are located between 6.4 keV and 6.9 keV , ranging from neutral 
iron to helium or hydrogen-like iron. For both sources these lines are expected at energies 
around 2 keV in the observer's frame, where the Chandra HETGS has a large collecting 
area and high resolving power. No significant line feature was detected in the spectrum of 
either source at a 3a level. In table 3 we list the upper limit equivalent width (EW) at two 
energies, 6.4 keV and 6.9 keV , for both MEG and HEG spectra for two assumed intrinsic 
line widths, (7 = and cr = 0.1 keV , respectively. The instrumental line response function 
is taken into account. Since the Chandra HETGS has better resolving power than ASCA, 
we can put a much tighter constraint on the equivalent width of the iron emission lines. We 
should emphasize that the low value of cr = is adopted to calculate the theoretically low 
limit of detectable EW and since we have no constraint on the upper limit of the intrinsic 
line width, the EW limit could be worse by a factor of 2 if cr = 0.5 keV . 

Y99 reported that they found a significant emission feature around 5 keV in the ASCA 
spectrum of PKS 2149-306. They proposed that this line was blue-shifted Fe-K emission 
with a bulk-motion velocity around 0.73c (head-on) and a rest-frame EW of 300 ± 200 eV 
(assuming an intrinsic line width of 0.01 keV , errors are 90% confidence). We do not detect 
any emission feature around 5 keV at 3cr level in either HEG or MEG spectra, and set an 
upper limit of the EW as low as 67 eV (Table 3). If we adopt the value of cr = 0.01 keV as 
used by Y99, the rest-frame EW hmit is 75 eV and still much smaller than lOOeF, which is 
their 90% confidence lower bound. A recent observation by BeppoSAX (Elvis et al. 2000) 
also addressed the existence of such a feature. 
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3.2.2. Absorption Lines 

X-ray resonant absorption lines are produced by highly ionized heavy elements. The 
detection of any absorption line depends on three factors: (1) the spectral intensity of the 
background quasar; (2) the physical properties of the intervening absorbers, specifically, the 
ion column density and velocity dispersion (both turbulent and thermal broadening); (3) the 
instrumental properties, such as the resolving power and effective area. Given that (1) and 
(3) are already known, we can constrain the properties of the absorbers. For demonstration 
we select the the spectrum of S5 0836+710 and calculate the expected signal-to-noise ratio 
(S/N) for Si XIV and O VIII (see figure 5), given the column density of Si XIV and O VIII 
at different redshifts and velocity dispersions {b, in units of km s~^). Two values of b- 
parameter, 200 and 500 km s~^, are adopted to bracket the typical velocity dispersions of 
hot gas associated with groups and clusters of galaxies. Since the strongest resonance lines 
of Si XIV and O VIII, the Lya transitions, are located at 2.01 keV and 0.654 keV (rest- 
frame), respectively, we calculate the S/N of Si XIV based on HEG (HEG has a better energy 
resolution there) and the S/N of O VIII based on MEG (HEG has no effective area below 
0.8 keV). The maximum redshift of O VIII is 0.3 because MEG has neghgible effective area 
below 0.5 keV. The horizontal line in figure 5 is the 3a detection level. The column density 
of O VIII or Si XIV must be at least 8 x 10^^ cm~^ or 5 x 10^^ cm~^, respectively, to produce 
an absorption line at or above the 3a level. 



4. X-ray Constraints on the Intergalactic Medium 
4.1. X-ray Forest 

The detection of the X-ray resonant absorption lines is closely related to the overall 
baryon density (Qb) and matter density (Qq) in the Universe. For a given level of metal 
enrichment, higher baryon density in general implies larger amounts of highly ionized metals 
and a higher possibility of detecting the X-ray absorption fines created by these metals 
in the spectrum of a distant quasar. On the other hand, as we predicted in Fang, Bryan 
& Canizares (2000), most of the high column density ions (specifically, at column density 
higher than 10^^ cm~^) are located in the collapsed, virialized objects, such as groups or 
clusters of galaxies, in the form of hot intracluster or intragroup gas^. Since flo governs 



^As we will discuss in the next section, numerical simulations show that the X-ray absorption material 
may also arise from filaments connecting the virialized structures; however, the resonant absorption lines 
produced by the low-column density ions (< 10^^ cm~^) from these filaments are generally too weak to be 
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the overall amplitude of the density fluctuations, changing Qq will dramatically change the 
spatial density of virialized objects, and so the distribution of ions at high column density. 

To investigate how Qj, and Qq can affect the X-ray spectra of high redshift quasars, we 
calculate the X-ray forest distribution function (XFDF, for details, see Hellsten, Gnedin & 
Miralda-Escude 1998; Perna & Loeb 1998; FC). The XFDF, d'^P/dN'dz, is deflned as the 
number of absorption systems along the line of sight with a column density between and 
N^ + dN^ per unit redshift, where iV* = iV(X*) is the column density for ion X*. For a source 
at redshift zo, the total absorption line number for ion X*, with a column density above a 
certain number Nq is given by 



Since the XFDF is determined by cosmological parameters, such as Qq ^-nd Qb (see figures 
2,3,4 of FC), by comparing the expected absorption line number (n) with observations, we 
may set the constraints on these cosmological parameters. 

Based on the Chandra observations we set the upper limit of different ion column den- 
sities (see section 3.2.2), which gives Nq in equation (1). The minimum detectable column 
density for O VIII is A^^ = 8 x 10^^ cm^^ and for Si XIV is A^^ = 5 x 10^^ cm^^. We find, 
even in the extreme case, where D,q — 1 and Zq — 1, the expected absorption line number 
for O VIII absorption lines is only n = 0.65 and even less for Si XIV. Here Zq is the metal 
abundance relative to the solar value where Zq = 1 indicates solar abundance, and wc inte- 
grate equation (1) to = 0.3 for O VIII and to = 2.17 for Si XIV. The Poisson statistic 
shows that the probability of detecting less than one O VIII absorption line is 0.53 given 
such an expectation. So with the current observations we are not able to constrain those 
cosmological parameters. This is also evidenced in FC. For example, in figure 7 of FC, even 
for the SCDM model, the expected number of O VIII absorption fines with N'^ > 10^^ cm~^ 
and z < 0.3 is equal to or below one. 

This method illustrates the possibility of constraining the cosmological parameters by 
observing more X-ray bright quasars with longer exposure time. For example, if a quasar is 
brighter than S5 0836-1-710 by a factor of two, by increasing the exposure time to 100 ksec, 
the minimum detectable column density for O VIII can be as low as A^q = 3 x 10^^ cm~^. For 
a low density cold dark matter model (LCDM) with Qq — 0.3 and = 0.7 ^, the expected 
number of O VIII absorption fines is 0.42, given a reasonable cluster gas fraction (~ 0.2) 



detectable with Chandra and even Constellation-X. 

^Here Oa = ^/Hq where A is the cosmological constant 




(1) 
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and metallicity (~ 0.5Zq for oxygen). The Poisson statistic shows that the probabihty of 
observing at least one hne is 95% if the expectation is three. This means by observing ~ 7 
sources which are twice brighter than S5 0836+710, we should have 95% chance of detecting 
at least one O VIII absorption line with an exposure time of 100 ksec. Previous observations 
show that at low and moderate rcdshifts, there are several tens of quasars having a flux 
higher than 5 x 10^^^ ergs cm^^s^^ between 2 — 10 keV, or 2 x 10'^^ ergs cm^^s^^ between 
0.1 — 2.4 keV (see, e.g., Brinkmann, Yuan & Siebert 1997; Yuan et al. 1998; Reeves & 
Turner 2000) . By increasing the number observations of bright quasars at low and moderate 
redshifts, there is a significant chance of detecting oxygen absorption lines. 

Similar constraints can be made with current and future missions, such as XMM and 
Constellation-X. For example, XMM and Constellation-X are able to detect O VIII at a 
column density as low as 7.9 x 10^^ cm"^ and 10^^ cm"^, respectively (see PC). Given such 
a low column density, we will detect as many as eight or ten O VIII absorption lines for 
a moderately luminous quasar at redshift of 2; ~ 1, respectively. Eventually this method 
may help us to understand the distribution of the baryonic matter in the highly overdense 
regions. 

4.2. X-ray Gunn-Peterson Test 

Recently large-scale cosmological hydrodynamic simulations have shown that a large 
amount of the hot/warm IGM (70% ~ 80%) lies in the regions at overdensity 5 < S < 200 
and temperature 10^ < T < lO''' (see, e.g.. Gen & Ostriker 1999; Dave et al. 2000). Here 
S = p/ (p) — 1 and (p) is the mean matter density. These regions are typically the extended 
filamentary structures connecting regions of high mass concentration. The X-ray Gunn- 
Peterson test may be applied to examine the possible "absorption trough" in the X-ray 
continuum spectrum of a distant quasar, produced by these diffuse structures. Shapiro & 
Bahcall (1980) and Aldcroft et al. (1994) discussed the absorption produced by a uniformly 
distributed IGM. Markevitch (1999) proposed that the X-ray Gunn-Peterson test can be used 
on a particular giant filament along the line of sight, in the direction of Aquarius. Based 
on numerical simulations and observations of the Lya forest systems wc know that instead 
of a uniform distribution, these filaments are actually clumpy structures in the universe. 
However, to simplify, we assume a uniformly distributed IGM and constrain its properties 
following Shapiro & BahcaU (1980) and Aldcroft et al. (1994). 

Gonsider the X-ray spectrum of a background quasar located at redshift Zq, the frequency 
of a photon emitted at vq will he 1/ — z/o(l + when the photon reaches the Earth. The 
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optical depth at frequency v (observer frame) is (Spitzer 1978) 

Z"^" dl 



(2) 







where ni{z) is the comoving number density of absorbing ion X* at a redshift of z and Ui is 
the ion absorption cross section at frequency v' — v z). The path length is dl. 

There are three dominant absorption processes contributing to the optical depth in 
equation (2): resonant absorption due to the atomic transitions, continuum absorption (edge) 
due to the bound-free transitions, and electron scattering. Electron scattering is independent 
of energy in the X-ray band and the only effect is uniformly lowering the X-ray continuum 
without altering the spectral shape. We ignore this effect because the absorption produced 
by this process will not be detectable without knowledge of the intrinsic quasar spectrum. 

To simplify we assume that the IGM is uniformly distributed hot gas in coUisional 
equilibrium, at temperature T > 10^ K. In such an environment, the hydrogen- and helium- 
like O, Ne, Mg, Si, and Fe XVII- Fe XXVI arc the most abundant ions (for ionization fractions 
in collisional equilibrium, see Mazzotta ct al. 1998). For the continuum absorption we include 
all the K-cdgcs (i.e., 18 K-cdgcs) and L-cdgcs for all the ions with L-shell electrons (i.e., 13 L- 
cdgcs). Wc obtain the threshold energies and cross sections from Vcrncr & Yakovlcv (1995). 
For resonant absorption, we select those transitions with oscillator strength greater than 0.1. 
This results in a total of 65 hues. The energy level and oscillator strength for each resonance 
hne are from Verner, Verner & Ferland (1996). The total optical depth Ti, then can be 
calculated based on equation (2). 

The quasar spectrum fy can be described as (Aldcroft et al. 1994) 



Here A is a normalization constant and / {Nh) is the the gas absorption where Nh is the 
hydrogen column density. The exponential term is the absorption of the IGM at a baryon 
fraction Vti, and temperature T. We first fit the observed spectrum with ilf, = 0. Then for 
each VLf) and T we fit the observed spectrum with equation (3) by varying A, Nh and F. We 
then calculate the statistic 



The Ax^(06,T') should follow the c/ii-square distribution with two degrees of freedom. In 
this way we can set confidence levels on Q;, and T. 




(3) 



(4) 



We select the MEG spectrum of S5 0836-1-710 because the MEG covers the entire range 
from oxygen to iron. Figure 6a displays the 3(7 contour of Ax^(Q6, T) for a flat, low-density 
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cosmological model (LCDM, Qq = 0.3, Qa = 0.7). The thick hnes indict an abundance 
(Zq = 0.5) and the thin hnes are for a low-abundance IGM {Zq — 0.1). For comparison 
we also show the constraint from the distortion of the cosmic microwave background (CMB) 
spectrum (which constrains the Compton y-parameter, the dotted line) as measured by the 
COBE satellite (Mather et al. 1994; Fixsen et al. 1996). The regions to the left of the sohd 
and dotted lines are allowed regions. 

We also display the total baryon density deduced from standard big bang nucleosynthesis 
(BBN) and the cosmic microwave background (CMB). BBN imphes n^h'^Q = 0.0388 ±0.0037 
(95% confidence level, hereafter cl) (see, e.g.. Buries, Nollctt & Turner 2000). The re- 
cent measurements of CMB anisotropy by BOOMERANG and MAXIMA gave Qbhjo = 
0.065lo:o}6 (95% cl) (Jaffe et al. 2000). We notice that while the fib measured in this way 
reflects the small angle (subdegree) anisotropies, the constraints based on the the Compton 
y-parameter reflects the CMB anisotropies at larger scales. We display both BBN and CMB 
results and the shadowed area shows the 95% cl region. 

From the figure we can find that while for a low-abundance IGM the X-ray Gunn- 
Peterson test is not able to put strong constraints on the properties of the IGM, for a 
uniformly enriched IGM with Zq — 0.5, marginal constraints can be made by combining 

the the Gunn-Peterson test and the Compton y-paramctcr from the CMB anisotropics. For 
example, for the CMB predicted baryon density, the temperatures of the IGM are constrained 
in a range with 10^'^ < 10^'^ K: temperatures above this range arc ruled out by the CMB 
spectrum, and at a temperature below this range the IGM will produce strong distortion in 
the quasar spectrum. In figure 6b we illustrate such a distortion in the MEG spectrum of 
S5 0836-1-710, based on a LCDM model. The dotted fine represents the absorbed spectrum 
by an IGM with Zq = 0.5 at T = 10^ K and Q^/i^o = 0.08. Given such an IGM, the model 
severely deviates from the observed spectrum. 

Based on the ROSAT PSPC observations of three high-redshift quasars, Aldcroft ct al. 
(1994) concluded that the X-ray opacity of a hot diffuse IGM is too small to constrain the 
IGM. For instance, they found that their X-ray Gunn-Peterson test did not constrain the IGM 
temperature above 10^ K with Qj, — 0.06 and solar abundances (Figure 11 in their paper). 
Our analysis indicates that the constraint can be improved significantly even with half-solar 
abundances, due to the high energy resolving power of Chandra HETGS. Our results are 
similar to their predictions with Chandra Low-Energy Transmission Grating Spectrometer 
(LETG), based on a simulated observation of 3C 273. However, due to the limitation of 
the small effective area, we are not able to constrain the IGM with low metallicity based on 
current data. In future, Constellation-X will increase the effective area by two magnitudes 
over Chandra HETGS and may eventually reveal the properties of the IGM. 
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5. Summciry 

In this paper we report the Chandra observations of two high-redshift, radio-loud 
quasars, S5 0836+710 and PKS 2149-306. 

1. The photon index (F) for S5 0836+710 is 1.39, consistent with previous observations 
with ROSAT and ASCA. Compared with previous observations, PKS 2149-306 has a 
rather soft spectrum with V — 1.255. Both photon indexes are consistent with the fact 
that RLQs have flatter spectra compared with RQQs. We also find excess continuum 
absorption (above the Galactic value) in S5 0836+710. The flux of S5 0836+710 is 
higher than in previous ASCA observations and puts this source as one of the most 
luminous objects at high redshift. 

2. No significant absorption or emission feature is found in both sources at or above ±3(7 
level. We put constraints on the possible emitters or absorbers. For example, Fe K 
emission hues were found in several RLQs (Reeves & Turner 2000), although not as 
frequently as in RQQs. We put tight upper limits on EW, as low as ~ 10 eV in both 
sources. We do not find the emission feature reported by Y99 in the spectrum of PKS 
2149-306 around 5 keV. We also examine the possible intervening systems and give 
upper limits on the absorbers' ion column density. 

3. We discuss the possibility of constraining the properties of the IGM. We develop a 
method to constrain cosmological parameters based on the X-ray forest method. How- 
ever, we are not able to set any meaningful constraints based on current data. The 
X-ray Gunn-Peterson test was also performed and we find that, by combining with 
the constraints from the CMB Compton y-parameter, this method can be applied to 
constrain a uniform, enriched IGM. 
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Fig. 1. — Isophotal contour map for the zeroth-order images of S5 0836+710 (left) and PKS 
2149-306 (right). The region size is 50" x 50". The images arc smoothed with a Gaussian 
{a = 1.2") for visual clarity. Both images can be well fitted by a single Gaussian with 
a ^ 0.38" for S5 0836+710 and a ^ 0.35" for PKS 2149-306. 

Fig. 2.— HEG spectra of S5 0836+710 (top panel) and PKS 2149-306 (bottom panel). The 
solid line and dotted line in each panel are the observed and fitted spectrum, respectively. 
The spectral model is an absorbed power- law and the best- fit parameters are from Table 2. 

Fig. 3.— MEG spectra of S5 0836+710 (top two panels) and PKS 2149-306 (bottom two 
panels) in wavelength space. The bin size is 0.041 for S5 0836+710 and 0.06A for PKS 2149- 
306. The observed spectra are first fitted with an absorbed power-law (same parameters as 
in HEG spectra), the residual is then fitted with a 5-order polynomial. The bottom figure 
in each panel shows the x, and the two dotted lines indicate the ±3(7 level. 

Fig. 4. — The cumulative distribution of the Poisson Probabilities fpr S5 0836+710 (solid 
line), PKS 2149-306 (dashed line) and the uniform distribution (dotted fine). The KS test 
shows that the null hypothesis that they are the same distributions is accepted at a significant 
level of 0.59 for S5 0836+710 and 0.32 for PKS 2149-306. 

Fig. 5. — Given the spectrum of S5 0836+710, we calculate the signal-to-noise ratio of Si XIV 
and O VIII absorption lines with different redshift {z) and velocity dispersion (6, in units 
of km s~^). Si XIV is based on HEG spectrum and O VIII is based on MEG spectrum. 
Since MEG has no effective area below 0.5 keV, the maximum redshift of O VIII is 0.3. 
The horizontal line is the 3(7 level. To obtain a Si XIV or O VIII absorption line in the 
spectrum of S5 0836+710 at 3a level, the ion column density must be at least 5 x 10^^ cm~^ 
and 8 x 10^^ cm~^, respectively. 

Fig. 6. — (a) Constraints on Qj, and the IGM temperature from the X-ray Gunn-Peterson 
test, based on the MEG spectrum of S5 0836+710. The top panel shows the 3a contour of 
(Qft/t^o) for each IGM temperature. The thick line indicts an enriched IGM with Zq — 0.5 and 
the thin line is for a low-abundance IGM {Zq = 0.1). The dotted line shows the constraint 
form the distortion of the CMB spectrum (the Compton ?/-parametcr) as measured by the 
COBE. The regions to the left of the solid (LCDM) and dotted (COBE) lines are acceptable 
(indicated by arrows), (b) The solid line is the observed MEG spectrum of S5 0836+710. 
The dotted line is the distorted spectrum absorbed by an enriched IGM {Zq = 0.5) with a 
temperature of 10^ K and fift/i^Q = 0.08. 
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Table 1. Chandra Observation Log 



Observation ID 


Object 


Redshift 




Start Date 


Exposure Time {ksec] 


01450 


S5 0836+710 


2.170 


2.83 


10/17/99 


58.7 


00336 


PKS 2149-306 


2.345 


1.91 


11/18/99 


37.3 


01481 


PKS 2149-306 


2.345 


1.91 


11/20/99 


55.5 



'^Galactic absorption from Murphy et al. (1996) (S5 0836-1-7104) and Dickey & Lockman 
(1990) (PKS 2149-306), in units of 10^° cm-^. 



Table 2. Chandra Spectral Fit 



Object 


Nh (1020 cm-2) 


Photon Index (F) 
(0.5 - 8 keV) 






S5 0836+710 
PKS 2149-306 


7.0 + 1.2^^ 
1.9 


1.388 + 0.012 
1.255 + 0.020 


3.96lE5:r, 


1.072/605 26.42 47.43/1^ 
0.965/605 7.58 46Mh^^ 



"^Normalization at 1 keV (observer's frame) in units of 10~^ photons cm~^s~^ keV~^ 
^Absorbed flux between 2 and 10 keV (observer's frame) in units of 10"^^ ^j-gg cm~'^s~^ 
'^Intrinsic luminosity between 2 and 10 keV (quasar frame) in units of ergs s~^, qo = 0.5 
^AU errors are quoted at 90% confidence 
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Table 3. The Maximum Equivalent Width of Possible Lines ^ 







HEG 




MEG 




Object 


Energy*^ a — 0'^ a — 


0.1 


(7 = (7 = 


: 0.1 


S5 0836+710 


6.4 


9 


72 


8 


49 




6.9 


14 


100 


15 


84 


PKS 2149-306 


6.4 


10 


83 


8 


53 




6.9 


12 


94 


12 


69 




17.0 


67 


188 


97 


215 



'^The equivalent width is measured in the rest frame of quasar 
(eV) 

''Line energy, in quasar frame (keV). The Fe K emission lines 
are located between 6.4 keV and 6.9 keV, the 17.0 keV line is 
the redshifted Fe-K line reported by Y99. 

■^Intrinsic line width (keV) 



S5 0836+710 



PKS 2149-306 




Right ascension 



Right ascension 
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> 


0.0100 


CD 












''cm' 


0.0010 




















0.0001 



S5 0836+710 




0.8 0.9 1.0 



2.0 

Energy (keV) 



5.0 



8.0 



0.0010 



X 0.0001 



PKS 2149-306 




0.8 0.9 1.0 



2.0 

Energy (keV) 



5.0 



8.0 
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- z = 0.0, b = 200, OVIII 
z = 0.3, b = 200, OVIII 

- - z = 0.0, b = 500, OVIII 

- z = 0.3, b = 500, OVIII 

- z = 0, b = 200, SiXIV 
z = l,b = 200, SiXIV 

- z = 0, b = 500, SiXIV 

- z = l,b = 500, SiXIV 



/ 



/ / 

/ 

^ ' ' ^^^^^^^^^ 



10 



15 



10 



16 



10 



17 



N (cm"') 



10 



18 




0.6 0.81.0 2.0 5.0 7.0 

Energy (keV) 



